Competitions in layered ruthenates: ferro- vs. antiferromagnetism and triplet vs. 

singlet pairing 
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Ru based perovskites demonstrate an amazing richness in their magnetic properties, including 
3D and quasi-2D ferromagnetism, antiferromagnetism, and unconventional superconductivity. Ten- 
dency to ferromagnetism, stemming from the unusually large involvement of O in magnetism in 
ruthenates, leads to ferromagnetic spin fluctuations in Sr2Ru04 and eventually to p-wave super- 
conductivity. A related compound Ca2Ru04 was measured to be antiferromagnetic, suggesting a 
possibility of antiferromagnetic fluctuations in Sr2Ru04 as well. Here we report first principles cal- 
culations that demonstrate that in both compounds the ferro- and antiferromagnetic fluctuations 
coexist, leading to an actual instability in Ca2RuC>4 and to a close competition between p-wave and 
d-wave superconducting symmetries in Sr2Ru04. The antiferromagnetism in this system appears to 
be mostly related with the nesting, which is the strongest at Q ~ (27r/3a, 2-7r/3a, 0). Surprisingly, 
for the Fermiology of Sr2Ru04 the p-wave state wins over the d-wave one everywhere except in close 
vicinity of the antiferromagnetic instability. The most stable state within the d-wave channel has 
vanishing order parameter at one out of three Fermi surfaces in S^RuCU, while in the p channel its 
amplitude is comparable at all three of them. 



In the last few years an understanding emerged, thanks 
to the progress in the high-T c problem, that s-wave and 
d-wave pairing are not entirely antagonistic. It is possi- 
ble to obtain either symmetry in the framework of one 
and the same model, depending on the actual values of 
parameters. It is, however, commonly believed that the 
triplet (p) and singlet (s,d) pairings are so different in 
their nature that more than just changing numerical pa- 
rameters of a model is needed to switch between these 
symmetries. Here we show that a realistic model of the 
superconducting layered ruthenate Sr2Ru04 is in fact un- 
stable with respect to both p- and d-wave pairing and it 
is a close competition between the two that determines 
the actual ground state. 

It was suggested a few years ago that Si - 2Ru04 may 
be a p-wave superconductor |lj . The main consideration 
was that the sister 3D compound, SrRuC>3, is a strong 
ferromagnet, so one could expect substantial ferromag- 
netic spin fluctuations in Sr2Ru04. It was also known 
that superconductivity in the canonical triplet supercon- 
ductor, 3 He, is due to ferromagnetic spin fluctuations, so 
it was natural to conjecture that the pairing in Sr2Ru04 
was also triplet. At that time there was hardly any ex- 
perimental evidence and no microscopic calculations to 
support this idea. Since then convincing experimental 
evidence has been collected (for review, see Ref. 0]) that 
the superconductivity in S^RuCU is indeed unconven- 
tional (not s-wave) and most likely triplet. Microscopic 
calculations revealed the mechanism for ferromagnetism 
in SrRuC>3 J3|,[| and demonstrated that a tendency to 
ferromagnetism is still present in Sr2Ru04, although it 
is weaker and does not result in an actual magnetic in- 
stability ||. 

The recent discovery of antiferromagnetism in 
Ca2RuC>4 forces us to reconsider this simple picture. This 



compound differs from S^RuCU only in that the RuOg 
octahedra are tilted and rotated [pi, as it is common in 
perovskites (cf. I^CuCU). This causes some modifica- 
tion of the hopping amplitudes compared to Sr2Ru04, 
but this modification has a tremendous effect on mag- 
netic properties: the material becomes antiferromagnetic 
with a substantial (> I ^/Ru) magnetic moment. This 
suggests that there should be some latent tendency to an- 
tiferromagnetism in Sr2RuC>4 itself, and casts doubt on 
the basic assumptions of Refs. |l],[| and others that the 
spin fluctuations in this material are predominantly of 
the ferromagnetic type. To answer this question one can- 
not rely upon analogies with other materials, but needs 
quantitative (at least semiquantitative) calculations. 

Fortunately, and unlike cuprates and most 3d-oxides, 
the conventional local density approximation (LDA) pro- 
vides a very good description of magnetic properties of 
the ruthenate-based perovskites. Previously we have per- 
formed calculations [||j] for various ruthenates for which 
crystal structure and magnetic properties are known ex- 
perimentally, and we found excellent agreement with 
the experiment: SrRuC>3 comes out ferromagnetic with 
the total magnetization 1.59 /zs/f.u. (experiment: 1.6), 
CaRuC>3 is a paramagnet on a verge of ferromagnetism, 
and the double perovskite Sr^RuYOe is antiferromag- 
netic with 3 /is/f.u., again in accord with the experi- 
ment. Finally, in Sr2Ru04 the paramagnetic state comes 
out more stable in the LDA calculations than either ferro- 
or antiferromagnetic one. 

The reason for ferromagnetism in SrRuC>3 (and near- 
ferromagnetism in CaRuOs) is now well understood [0: 
There is substantial oxygen density of states at the 
Fermi level in these ruthenates (due to strong p — d 
hybridization), and the difference between the FM and 
the AFM state is that in the latter case the oxygen 



is not spin-polarized. The oxygen ion has considerable 
Stoner energy, which is lost in the AFM case. Fur- 
thermore, this additional energy is entirely lost for the 
q = {7r,7r,7r} AFM ordering, two thirds of that is lost 
for q = {7r,7r,0}, and one third for q = {it, 0,0}, com- 
pared with the FM ordering (q = {0,0,0}). This al- 
lows one to construct a g-dependent Stoner interaction, 
I(q) ~ 0.46 eV/(l + 0.08q 2 ), where q is measured in units 
of n/a. This interactions strongly favors a FM instability, 
and whether or not the actual instability occurs depends 
on the density of states at the Fermi level, according to 
the Stoner criterion, I(0)N(0) > 1. It appears that in 
SrRuOa this condition is satisfied, IN — 1.23, and the 
material is a FM. In CaRuOs the smaller ionic radius of 
Ca leads to a smaller Ru-Ru distance and thus to larger 
distortion. A peak in the density of states that exists in 
SrRuC>3 is washed out and the material is on the border 
line, 77V » 1. 

The same mechanism is operative in S^RuC^ : for an 
individual RuC>2 plane we obtain a (2D) Stoner factor 
I(q) ~ 0.43 eV/(l + 0.08q 2 )_, favoring ferromagnetic spin 
fluctuations in the plane ||. However, the 2D charac- 
ter of the band structure of Sr2Ru04 introduces addi- 
tional complications. As discussed in Refs. 0,H|, of the 
three Fermi surface sheets one (7) is quasi-isotropic 2D, 
and two (a and /3) are quasi-lD. The latter can be vi- 
sualized (cf. Fig.l in Ref. M) as a system of parallel 
planes separated by Q — 2n/3a, running both in the 
x and y directions. This is true in the nearest neigh- 
bor ddn tight-binding model, while in reality due to the 
next hoppings the planes are warped and reconnected 
at the crossing lines, to form two pseudo-square prisms, 
obtained in LDA calculations and observed experimen- 
tally Naturally, such a Fermi surface should give rise to 
sizable nesting effects at the wave vectors k = (Q,k y ), 
k = (k x ,Q), and especially at k = Q = (Q,Q). This 
would lead to AFM spin fluctuations at these vectors, 
in addition to the FM fluctuations discussed above. To 
check, we have integrated the LDA band structure of 
Si'2Ru04 to get the bare RPA susceptibility, 
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where / is the Fermi distribution function, i and j label 
the three bands. All |kz > states were classified accord- 
ing to the maximal ti g character, xy, yz, and zx, and the 
matrix element M is taken to be 1 between two states 
which have the same maximal character and otherwise. 
This is, of course, a rather crude approximation, but it 
should reveal the qualitative behavior of xa- The results 
are shown in Fig.l. Roughly speaking, 



Xo(q)=7V(0)+x«(q), 



(2) 



1 - J(«)xo(q) 1 - i(q)N(o) - J(«)x n (q) ' 



(3) 



This form implies two different kinds of spin fluctua- 
tions: FM ones, q= 0, and AFM ones, at q = Q. If 
I(Q)N(0) +I(Q)Xn(Q) > 7(0)JV(0), the AFM fluctua- 
tions are stronger. This seems to be the case in S^RuCU : 
our calculations yield I(0)N(0) = 0.82, in good agree- 
ment with the experimentally observed susceptibility en- 
hancement, and I(Q)N{0) + I(Q)Xn(Q) = 1-02 (which 
actually corresponds to an instability with respect to 
tripling of the unit cell both in x and in y) . Since no insta- 
bility is observed in the experiment, nor in the direct cal- 
culations, we conclude that the approximate treatment of 
the matrix elements in EqlP leads to an overestimation of 
Xn by at the very least 2%, but the conclusion that AFM 
fluctuations are stronger or at least comparable with the 
FM ones likely holds. 
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where Xn is the nesting-dependent contribution. The to- 
tal susceptibility can then be expressed as 



FIG. 1. Calculated bare susceptibility for S^RuCU 

The direct way to test this experimentally is via neu- 
tron scattering B. There is, however, an indirect ar- 
gument in favor of strong AFM spin fluctuations. In- 
creasing the effective dimensionality by adding additional 
RuC>2 layers, one can increase N(0) and eventually get 
a FM instability. Experimentally this happens when the 
number of layers is 3 or maybe even 2 S . Another pos- 
sible (but not guaranteed) effect of adding layers is in- 
creased z-dispersion and thus deteriorated nesting. On 
the other hand, reducing the next-nearest-neighbor hop- 
ping should improve nesting and make an AFM transition 
more likely. One expects such a reduction from rotating 
the RuOe octahedra pi, as for instance in Ca2Ru04. 
Indeed, experimentally Ca2Ru04 is an AFM with a mag- 
netization of 1.2-1.3 /is/Ru and Tm ~ 150 K. Moreover, 



this AFM state is remarkably different from typical Mott- 
Hubbard insulators, driven by strong Coulomb correla- 
tions. First, although the conductivity grows with tem- 
perature, the functional dependence is consistent with 
a variable-range hopping and not with activation. Sec- 
ond, there is substantial density of states at the Fermi 
level, as evidenced by specific heat measurements. These 
two facts indicate that Ca2Ru04 is not a simple insula- 
tor, but a metal with disorder localized carriers (which 
is in turned helped by strong coupling between the spin 
and charge degrees of freedom B ) . We performed LDA 
calculations for Ca2Ru04 similar to those reported in 
Refs. ||]4|]and found a magnetic moment of rs 1.5 hb (of 
which « 1 fiB is inside the the Ru MT sphere and the 
rest mostly residing on the apical oxygens) for Ca2RuC>4 
(in agreement with the experimental M = 1.3 /is). In 
fact, we also find a FM instability, and that the FM and 
the AFM states Ca2Ru04 degenerate within the LDA 
to within a few meV/atom, indicating a close competi- 
tion between these two magnetic states pH| . For gen- 
uine Mott insulators the LDA either fails to reproduce 
the magnetic instability, or underestimates the magne- 
tization. The calculated density of states is sizable, 
N(Ep) = 1.6 st./eV spin f.u.), but the corresponding 
effective mass is large (the in-plane average |L2| n/m= 
7.5 x 10 20 cm~ 3 /TOo), because Ru is strongly spin polar- 
ized and thus Ru-Ru nearest neighbor hopping is sup- 
pressed with AFM order. This also facilitates localiza- 
tion. 
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a, weight of the nesting contribution 

FIG. 2. The maximum eigenvalue for the coupling ma- 
trix in the singlet and in the triplet channels, as a function 
of the relative strength of the AFM component in the bare 
electronic susceptibility. 



We now enumerate the relevant magnetic interac- 
tions in the (Sr,Ca) 2RUO4 system. The first two have 
been discussed above, and they are the Stoner ferro- 
magnetism (Q = 0), nesting-derived antiferromagnetism 
(Q « {2tt/3, 27r/3}). As usually, there is also superex- 
change (Q = {ir, ir}). One should be reminded at this 
point, that superexchange per se is not a strong correla- 



tion phenomenon. In a quasi-one-electron approach like 
LDA, for sufficiently large exchange splittings, an AFM 
state is lower in energy than the FM one by roughly 
Zt 2 /A, where Z is coordination, t is the hopping, and 
A is the exchange splitting. Since t is usually overesti- 
mated in LDA, and A is usually smaller than Hubbard 
U, the effect of superexchange is usually overestimated, 
not underestimated in LDA (see Ref. M for more). 

While superexchange is a universal mechanism and 
should also be operative in Ca2Ru04, it hardly plays a 
leading role. Otherwise, the LDA calculations that gave 
the right Ru moment would have overstabilized the AFM 
state, while in reality the opposite is true (they come out 
degenerate instead of the AFM being the ground state). 
The most likely cause for the instability is nesting. Fur- 
thermore, superexchange is of less relevance for Sr2Ru04 
because it manifests itself for finite amplitude spin fluc- 
tuations, and not in the low-energy spectrum. 

Assuming that the LDA gives a reasonable description 
of the spectrum of spin fluctuations in Sr2Ru04, we can 
discuss the consequences that the AFM spin fluctuations 
may have on superconductivity. Recall that in p-wave 
superconductivity only the small q spin fluctuations are 
pairing, while in the <i-wave case they are mostly pair- 
breaking. For a spectrum with a complicated q depen- 
dence, as the one described by Eq||, the most favorable 
state is defined by an interplay of the Fermi surface ge- 
ometry and the structure of the effective interaction. In 
the weak coupling near T c the gap equation looks like 



Au 



k'j 



I^ki,k'jAk' 



(4) 



where for the singlet (d) pairing A is the order param- 
eter, and Vki,k'j is (negative) pairing interaction, which 
in the same approximation as the one used in Ref. H is 
given by 



V(q = k - k') = - 



/ 2 (g)xo(q) 

1-P(?)x8(q)' 



(5) 



For the triplet (p) pairing A is the amplitude of the order 
parameter, usually denoted as d, and Vki,k'j now includes 
the sign-changing angular factor: 



V(q = k - k') = 



v k -v k / / 2 (g)xo(q) 



(6) 



«kWk' i-/ 2 0?)xo(q) 

The largest eigenvalue of the matrix V in Eq.[| defines 
the critical temperature, and the corresponding eigen- 
vector defines the anisotropy of the gap near T c . Cor- 
respondingly, the instability in the singlet channel is de- 
fined by the same matrix, but with the opposite sign (the 
interaction is repulsive in singlet channel), and without 
the angular factor Vk ' Vk ' . We have solved these eigen- 
value problems numerically using a discrete mesh of k 
points at the Fermi line. We used model susceptibil- 
ity (0), adding an adjustable reduction a for the nesting 



part, Xo(q) = N(0) + ax„(q), and look for the solutions 
at different a < 0.98 (i.e., below actual AFM instabil- 
ity). Furthermore, to simplify the calculations, we used 
an analytical form for \ n (q) that recovers the main qual- 
itative characteristics of the numerical result of Fig.l, 
namely x„(q) = A[cos(aq x ) + cos(aq y )] + B[cos(2aq x ) + 
cos(aq y )] + C[cos(aq x ) + cos(aq x )]. 

The results are shown in Fig. 2. Amazingly, only 
in close vicinity to the actual AFM instability (.85 < 
a < 0.98) the superconducting transition happens in 
the d-wave channel. For smaller a the p-wave state is 
favored. An interesting question is, what is the angu- 
lar and/or interband anisotropy of the gap. We illus- 
trate this anisotropy for the critical value of a (Fig. 3). 
The p-wave state is relatively isotropic, indicating that at 
least in the Fermiology there is no reason for substantial 
interband anisotropy ("orbital-dependent superconduc- 
tivity"), suggested in Ref. [Q. The d-wave, to the con- 
trary, is substantially anisotropic beyond the standard 
anisotropy associated with the nodes along the {tt, tt} 
direction, namely the so-called a pocket of the Fermi 




FIG. 3. Relative magnitude of the order parameter in 
the singlet (left) and the triplet (right) states. Parameter 
a is the same for both panels (0.83), corresponding to the 
coupling constant in both channels A ~ 0.5. 

Our conclusions are as follows: (1) LDA calculations 
yield a self-consistent solution with the AFM ordering for 
Ca 2 Ru04, degenerate with a FM solution. The nonmag- 
netic state is considerable higher in energy, in agreement 
with the experiment. The calculated magnetization value 
of rs 1 fi B is in agreement with the experiment as well. 
(2) This fact suggests that besides the known tendency 
to ferromagnetism, there is an intrinsic tendency to an- 
tiferromagnetism in Ru based layered perovskites. (3) 
Analysis of the Fermi surface geometry of Sr2Ru04, as 
well as direct calculations, indicate a strong nesting at 
Q = {±27r/3, ±27t/3}, as well as a weaker nesting for the 
wave vectors connecting these four points (Fig. 1). This 
result also suggests that the AFM ordering in Ca2Ru04 
is mostly due to nesting, although superexchange may 
play some role as well. It is, however, unlikely that 
this compound is a Mott-Hubbard insulator. (4) The 
spin fluctuations in Sr2Ru04 have both a FM and an 
AFM component, of comparable magnitude, thus mak- 



ing d-w&ve superconductivity a strong competitor with 
the p-wave state. It is possible that the system may be 
driven to the d-wave superconductivity by an external 
force, for instance by pressure. (5) Spin fluctuation spec- 
trum combined with the actual fermiology of Sr2Ru04 
produced a p-wave state with little angular or interband 
anisotropy. (6) The sister compound Ca2Ru04 is a low 
carrier density metal, where carrier are localized due to 
disorder and electron-magnon scattering, and where an- 
tiferromagnetism is primarily due to the nesting effects. 

We acknowledge useful discussions with M. Sigrist, 
J.W. Lynn and E. Demler. Computations were per- 
formed using facilities of the DoD HPCMO ASC Center. 
Work at the Naval Research Laboratory is supported by 
the Office of the Naval Research. 
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Imai et al (preprint submitted to PRL) recently reported 
NMR relaxation rate measurements which they inter- 
preted as indication of the absence of AFM spin fluc- 
tuations. However, their data only relate with the spin 
susceptibility at q = (it, tt), while our calculations pre- 
dict spin fluctuations at q = Q =|(7T, tt). 
In the calculation (to be published elsewhere) the n — 2 
compound Sr3Ru207 is not yet ferromagnetic; experi- 
mental reports are contradictory, but the n = 3 material, 
Sr4Ru30io is positively FM. 

We do not discuss here the smearing of the density of 
states due to additional Bragg reflections in a distorted 
structure, an effect that controls the loss of magnetic or- 
dering from Sr_Ru03 to Ca.RuOs,. 

We used crystal and magnetic structure of Ca2Ru04 as 
recently measured by means of neutron scattering fef. 
The full-potential linear augmented plane wave method, 
as described in |M,k| was utilized. The ferromagnetic 
structure appeared lower than the AFM one by ~ 0.5 
mRy/Ru. While the technical accuracy of the computa- 
tions is better than that, the LDA itself is not as ac- 
curate in magnetic energy differences; LDA and GGA 
results, for instance, usually differ by a larger margin. It 
is not excluded that on top of the (tt, n) AFM structure 
reported in Ref. M there is an incommensurate super- 



structure with Q~ 0.6(7r/o, 7r/6). This would lower the calculations will be published elsewhere. 

energy slightly and probably render this AFM state more [12] The effective mass has substantial (factor of 3) in-plane 

stable than the FM one. A more detail account of these anisotropy. 



